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Abstract 

We study contributions of quark-squark-gluino interactions in Minimal Su- 
persymmetric Standard Model (MSSM) to B -> PP (PP = Kir, vrvr, K K ) 
decays using QCD improved factorization method for the evaluation of the 
hadronic matrix elements and taking into account renormalization group run- 
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ning of the Wilson coefficients from SUSY scale (~ rriq) to the low energy scale 
(= nib) applicable to B decays. Using the most recent experimental data we 
obtain constraints on flavor changing Supersymmetric (SUSY) parameters 
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(dij) ll,rr- For AS = —1 processes, b — > S7 is usually considered to give 
the strongest limits. We, however, find that in some part of the parameter 
space B — ► Ktt processes give stronger bounds. Implications for B® — Bg mix- 
ing is discussed. We also study constraints obtained from AS" = processes 
B — ► 7T7T, KK, /?7 and — B& mixing. In this case, in a large part of the 
parameter space B& — B& provides the best bound, but B — > K~~K°,pj can 
still give interesting constraints in a complementary region of the parameter 
space. 
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I. INTRODUCTION 



Recent results from 5-factories at CLEO, BaBar, Belle have attracted a lot of attentions. 
Data from these experiments provide new opportunities to further study the Standard Model 
(SM) and also to study physics beyond the SM. Rare charmless hadronic B decays, such as 
B — > PP(PP = Kir, -kit, KK), are particularly interesting. These processes being rare are 
sensitive to new physics beyond the SM. Minimal Supersymmetric Standard Model (MSSM), 
which is an extension of the SM emerges as one of the most promising candidates for new 
physics beyond the SM. MSSM resolve many of the potential problem of the SM, for example 
the hierarchy problem, unification of SU(3) x SU(2) x U{\) gauge couplings, and so on [1]. 

B — > PP decays have been studied extensively within the framework of naive factoriza- 
tion scheme. In the last few years, QCD improved factorization method have been developed. 
This method incorporates elements of the naive factorization approach (as its leading term) 
and perturbative QCD corrections (as subleading contributions) allowing one to compute 
systematic radiative corrections to the naive factorization for the hadronic B decays [2,3]. 
This method provides a better understanding of the strong interaction dynamics for B decay 
into two light mesons. Therefore results obtained using QCD improved method are expected 
to be more reliable compared with those obtained by using naive factorization. Lot of studies 
have been done in this direction so far [4,5]. 

In this paper we study implications of the recently measured B — > PP branching ratios on 
SUSY flavor changing parameters, specially arising from the quark-squark-gluino (q — q — g) 
interactions in the model [6]. Effects of SUSY interactions on various processes involving 
B mesons have been studied previously [6-11]. Here we will concentrate on the q — q — g 
interactions and re-analyze B — > PP decays with several improvements. Firstly we use the 
QCD improved factorization method to compute the two body hadronic B decays. Secondly, 
we improve the quark level effective Hamiltonian by including QCD running of different 
Wilson coefficients (WC) from the SUSY scale (= rriq) to the low energy scale (= m b ) which 
is applicable to B decays. Previous authors while doing similar analysis have neglected this 
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RG evolution of WC's. We find that this RG evolution of WC's have non-negligible effects 
on the branching ratio calculations. 

Using the most recent experimental data from different 5-factories on the B — > PP 
decays we constrain the flavor changing SUSY parameters arising from the quR—QjLR mixing 
due to the presence of q — q — g interaction vertices in different one loop contributions to 
B — > PP decays. Similar kind of squark mixing can also occur for other flavor changing B 
decays. For AS = — 1 processes, b — > 57 is usually considered to give the strongest limit 
on such mixing. We, however, find that in some part of the parameter space B — > Kix 
processes provide better bound. We also discuss the implications of B® — B® mixing in 
this context. Furthermore, we also study constraints on FCNC parameters obtained from 
AS = processes B — > tttc, KK, fry, and B b — B d mixing. We find that B° d — B® mixing 
provides the best limit in a large part of the parameter space, but in certain region of the 
parameter space B — > tttt, KK, p7 can still give better constraints. 

II. SQUARK-GLUINO CONTRIBUTIONS TO B — > PP DECAYS 

In the SUSY, squark and gluino induce large one loop contribution to B — > PP decay 
from penguin and box diagrams because they are related to strong couplings. We compute 
the effect of flavor changing contributions to B — > PP arising from q — q — g interaction 
vertices using the mass insertion approximation method [6,12]. In this method, the basis of 
the quark and squark states are chosen in such a way that their couplings with gluinos are 
flavor diagonal. The flavor changing current arises from the non-diagonality of the squark 
propagators. We denote the off-diagonal terms in the squark mass matrices (i.e. the mass 
terms relating squark of the same electric charge, but different flavor), by A, then we expand 
the squark propagators as series in terms of {5ij)AB = i^-ij) ab / 'm 2 where m is an average 
squark mass, A, B = L,R and i,j are the generation indices. The notation (5is)ab and 
($23) ab denote d — b and s — b mixing respectively. These two parameters signify the size of 
flavor changing interactions in SUSY through the q — q — g interactions. 
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At one loop level, for qi — qi mixing (LL type ), the new effective Hamiltonian relevant 
to B — > Kn, ixix decays arise from SUSY penguin and box diagrams with gluino-squark in 
the loops is given by 

KW = -^v^jzc^o,. (i) 

Here we have normalized the WC's according to the SM ones. The operator Oj is the same 
as SM definition with 

3 = qa^Lbi QjluWj, 04 = QiYLbj ^ q'^Lq'^ 

q< q> 

5 = qrfLbi X: q'jl.Rq'j, 6 = qtfLbj £ q'^Rq',, (2) 

q> q> 

where L(R) = (f ± 75), q — d, s, q' — u, d, s. The WC's are given by [6]: 

C " (M -> = H B ' W + 5 B2(I) + S Pl W + 5 P * W ) 

where j = 1,2 when q — d,s, respectively and x = rri~/m~ . The definition of the functions 
Bi(x), Pi(x) can be found in Ref [6]. 

The above WC's are obtained at the SUSY scale M susy . We take this scale as rriq. We 
then evolve these WC's from this high scale down to the scale m b using the RG equation. 
At the lower scale, after the RG evolution, the WC's take following form: 

C 3 LL (m*>) = b\C^ L {M susy ) + b 4 3 C^ L (M susy ) + b\C^{M susy ) + btC^ L (M susy ) 
C L , L {m b ) = blC^ L (M susy ) + b\Cl\M susy ) + b\C^\M susy ) + b\C^ L {M susy ) 
Ct L (m b ) = bl^ L {M susy ) + b\Ci L {M susy ) + blC^ L (M susy ) + btC^ L (M susy ) 
C^ L (m b ) = blC^ L (M susy ) + btC^ L (M susy ) + blCt L (M susy ) + ^C 6 LL (M SUSJ/ ) (4) 
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M susy (GeV) 


4 


bi 




4 


bl 


bi 


bl 


bl 


100 


1.144 


-0.222 


0.009 


0.086 


-0.322 


0.989 


-0.022 


-0.197 


500 


1.208 


-0.299 


0.018 


0.139 


-0.432 


1.002 


-0.042 


-0.298 


1000 


1.235 


-0.328 


0.022 


0.164 


-0.474 


1.009 


-0.051 


-0.342 


M susy (GeV) 


H 


H 


bl 


bl 


bl 


bi 




bl 


100 


0.011 


0.039 


0.925 


0.058 


-0.051 


-0.173 


0.330 


1.742 


500 


0.011 


0.050 


0.905 


0.083 


-0.064 


-0.258 


0.478 


2.084 


1000 


0.011 


0.053 


0.898 


0.094 


-0.068 


-0.294 


0.540 


2.227 



TABLE I. Values of different WC's after RG evolution for three values of M SUS y = 100, 500 and 
1000 GeV. 



bj are the constant after running, depend on the choice of the initial scale, as shown in Table. 
I. 

We find that the WC's with and without the RG evolution can be different by a factor 
of two in some cases. This implies that to have a reliable result and consistent computation 
one has to take into account the RG evolution of WC's. 

Apart from LL type of operators, one can also have operators with different chiral struc- 
tures, such as RR, LR and RL type. The operators and WC's for RR mixing can be obtained 
by the exchange L <-> R in Eqs. 2 and 3. LR and RL mixing cases can also be easily studied. 
However, in these two cases there are more stringent constraints from gluonic dipole moment 
operators for B — > PP decays. The details have been studied in [13]. We will not discuss 
these two cases further. 

III. DECAY AMPLITUDES FROM QCD IMPROVED FACTORIZATION 

To obtain the branching ratios of B — > PP decays, one needs to study the matrix 
elements, < PP\H e ff\B >. The usual approach is the naive factorization method, where 
the hadronic matrix element can be approximated as a product of two single current matrix 
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elements, then it is parametrized into meson decay constant and meson-meson transition 
form factor. 

Previous analysis of SUSY contribution to B — > PP decays have been performed using 
the above mentioned naive factorization method. However, in the present analysis we will use 
the QCD improved factorization approach [2] to compute the SUSY effects to the B — > PP 
decays. This formalism incorporates the basic elements of naive factorization approach as 
its leading term and perturbative QCD corrections as subleading contributions allowing one 
to calculate radiative corrections to the naive factorization for the hadronic two body decays 
of B meson. 

In our case the total effective Hamiltonian for B — > PP decays have two sources, the 
usual SM part H^ff and the SUSY part H^f . We use the SM results obtained in Ref. 
[2]. The results for the two types of SUSY contribution from LL and RR mixing can be 
easily obtained once one understands the SM contributions since the operators from SUSY 
interactions discussed here are all strong penguin type. In the LL mixing case, the total 
WC's are the sum of SM WC (C[^) and SUSY contributions. 

/-itotal ^SM | ^~<LL s-itotal ^SM , ^LL 

o 3 — o 3 -+- o 3 , <_/ 4 — (_/ 4 to 4 , 

Qtotal qSM _|_ qLL Qtotal qSM _|_ qLL 

For RR mixing, due to the different chirality of the new operators the total WC's in the 
final results are replaced by Cg^jg = C|^ 5 6 — C^f 5 6 . The two type of mixing, LL and RR 
have similar contributions. Knowing how one of them affects B — > Kn, nn decays, the other 
can be easily estimated. 

We write down the decay amplitude of B° — > it + it~, K~it + for an illustration, 

A{ B° -> vr+Tr-) = ^iMml - ml)F B ^(0) [K««i + < + <o 
+R*(a% + <)) + V cb V c * d {a c A + a c l0 + R n {a c 6 + a c 8 )) 

( "' ' ' 2 r K & K*A + (v ub v: d + v cb v* d ) (b 3 + 26 4 - hr + \bT) 



+l 7=2 fBf * 



(6) 



A(B° - ^-tt+) = -^Mm| - m 2 K )F B ^(m 2 K )[V ub V: s ( ai + < + < + + <)) 
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+V ch V;M + a c w + R K (a c 6 + a c 8 )} 

+^ifsfMv ub v: s + v cb v; s )(b 3 - l -bT) 



(7) 



Similarly for other B — > Kn, tvtt, KK modes. 

Here Rk = 2m 2 K /m s mb. cti and foj coefficients are related to the WC's. Including the 
lowest a s order corrections, af's are given by 

C2 Cf^ol s 



ai = C 1 + — 
a 2 = C 2 + — 



1 + — 7— ^Vk - ' 



I + 



47T 
47T 



+ 
+ 



iV c N c 



-H 



Kit 



N r N r 



H, 



ttK 



a\ = C, 



total I ^3 
4 



Qtotal 



V /-itotal 1 

a 7 = C 7 + — 



TV 

Qtotal 



N 
I + 



1 + 



1-6 



= C 8 + 

a 9 — Cg + 



P 

a 10 



C10 + 



CV 
iV 
C10 

c 9 



iV 



1-6- 
I + 

1 + 



An 



4tt 
(-K) 



+ 
+ 



Cfol 9 P p k,2 C\ otal C F na t 
Air N c N c N c 

Air N r 



H 



Kw 



C 8 C F ira s 
f — — - — (-n wK ) 



4tt 

C F Oi 



+ 



N c N c 

p p,EW 
9TT iV r 



4tt 



C10 Cfvtq; s 

+ ~7T7 77 J^nK 



47T 



A" 



+ 



N c N c 

p p,EW 

9tt N r 



C 9 C F na s 
+ — 77 — HKtt 



N c N c 



(') «■(') 



r is o, .£ ■ 



M 2 Mi5 J K,2> J K,Zi 1 K,2 



P. 



(8) 



,p,EW 



where C F = (iV c 2 - 1)/(2N C ), and iV c = 3. The quantities H 
and Pk^ W are hadronic parameters that contain all nonperturbative dynamics. The vertex 
corrections V M (M = ir, K) are given by 



V M = 12 In — - 18 + dxg{x)(t)M{x) 
fi Jo 

V' M = 12 In ^ - 6 + f 1 datf(l - x)0 M (:r) 
« J 



(9) 
(10) 



H 

g( x ) = 3(1—^- \n(x) - in) + [2L* 2 (:r) - ln 2 (a;) + - (3 + 2m) - (x <-> 1 - x)] . (11) 



1 -x 



1 -x 



where 4>m(x) is the leading- twist light cone distribution amplitude of light pseudoscalar 
meson [14]. This distribution amplitude can be expanded in Gegenbauer polynomials. We 
truncate this expansion at n — 2. 
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<j )M (x, f j,)=6x(l-x) l + a™(n)C{ 3/2) (2x-l) + a™{n)C ( 2 



,(3/2)(2x-l) 



(12) 



where C{ (u) = 3u and C 2 (u) = §(5w 2 — 1). The distribution amplitude parameters 
af 2 for M — K, n are: af = 0.3, a 2 = 0.1, c*f = and a 2 = 0.1. The detailed expressions 
can be found in [2]. 

For the penguin contributions, P Ki2 , ^a",3) Pk^ W an d Pk^ W are given by 



4 .7716 2 

- In (xK\Sr>) 

3 // 3 V pJ 



_|_ ^jtotal 



x 



i^J w = (d + iV c C 2 ) 

-Pk\3 = ^1 

P P kT = (Ci + iV c C 2 



1^/ In 1^ _ ( n/ _ 2 )Gx(0) - (a c ) - 



4 m 6 2 

3 ln 7 + r w 



4. m 6 2 - 



_|_ Qtotal 



-3C 12 fJE-Mx) 
Jo 1 — x 



i n 1^ _ (n/ _ 2 )G X (0) - G K ( Sc ) - G x (l) 



2C 



ii 



4, m 6 2 . 



3C 



12 



(13) 



with 



Gr(s) = / cLxGk(s — ie, 1 — x)(f)K{x) 
Jo 

Gk(s) — / cLxGk{s — ie, 1 — x)<f£(x) 
Jo 

Gr(s, x) = —4 / duu(l — u) ln[s — u(l — u)x] 
Jo 



(14) 



where <j> p (x) = 1. 



The hard spectator contributions to the coefficients Oj's are parametrized in terms of a 
single (complex) quantity H' MiM2 which suffers from large theoretical uncertainties related 
to the regularization of the divergent endpoint integral. Following [2] we use H nK = H KK = 
0.99 at the scale /i — m& and 



(15) 



All scale-dependent quantities for hard spectator contributions are evaluated at Hh = 
V^hfrib with A h = 0.5 GeV. 

Terms containing the coefficients 6j and bf are from annihilation contributions. The 
annihilation coefficients are 

&! = %C 2 A\, b 3 = ^L[Cl otal A\ + Cf al (Al + A{) + N c C l ° tal A{] 

■'"c c 

& r = SlO^i + C 7 {A\ + 4) + N C C S A{] 

c 

fcr = ^[Cio4 + C 8 ^]. (16) 

with 



^ I -2 V 2 



18 X A -4 + — \+2r z x Xj L 



4 = 127ra4(2Xl-X A ) (17) 

where r x = 2jj, h /m b , Xa = j^dy/y parameterizes the divergent endpoint integrals. We 
use Xa = ln(m B / A h ). All scale dependent quantities for annihilation contributions are 
evaluated at jih- 

In the QCD improved factorization method the weak annihilation contributions are power 
suppressed as AQ C D/ m b and hence do not appear in the factorization formula. Apart from 
this, they also show end-point singularities even at twist-two order in the light-cone expan- 
sion of the final state light mesons and therefore can not be calculated systematically in 
the context of hard scattering approach. One can bypass this problem, by treating these 
different end-point singularities as a phenomenological parameters. But this induces model 
dependence and numerical uncertainties in the calculation. 

Using the above technique we compute the different two-body branching ratios of B 
mesons within the framework of SUSY and compare the calculated results with the experi- 
mental data to constrain the different flavor changing parameters of SUSY. The experimental 
branching ratios used in this analysis are listed in Table II. We average the data from CLEO, 
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BaBar, Belle 5-factories as their results are uncorrelated. To impose the limit on these fla- 
vor changing parameters we require that the theoretically computed branching ratios should 
not deviate from the experimental measurement by 2a. If the experimental branching ratio 
has only upper limit, we select the most stringent one from the experimental data. It is 
clear from Table II that for B — > tt°tt and B — > KK modes have only the 90% C.L. upper 
limit. 



TABLE II. The branching ratios for B — >■ PP in units of 10 6 . 



Branching ratio and 


Cleo 


Belle 


Babar 


Averaged 


CP asymmetries 


[15] 


[16] 


[17] 


Value 


Br{B u tt-K°) 


IS Q+3. 7+2.1 
i8 -°-3.3-1.8 


22.0 ± 1.9 ±1.1 


17.5l{;^ ± 1.3 


19.7 ± 1.5 


Br{B u 7T°K-) 


19 q+2.4+1.2 
iz - y -2. 2-1.1 


12.8 ± lA±\i 


12.8±i;f ± 1.0 


12.8 ±1.1 


Br(B d -► TT+K-) 


io n +2.3+1.2 
i8 - U -2.1-0.9 


18.5 ± 1.0 ±0.7 


17.9 ±0.9 ±0.7 


18.2 ±0.8 


Br(B d -► 7r°K°) 


in O+4.0+1.7 
iz -°-3.3-1.4 


12.6 ±2.4 ±1.4 


10.4 ±1.5 ±0.8 


11.2 ±1.4 


Br(B u vr-vr ) 


a fi +1.8+0.6 
4 -°-1.6-0.7 


5.3 ±1.3 ±0.5 


5.51^ ±0.6 


5.3 ±0.8 


Br(S d vr+vr-) 


A r+1.4+0.5 
2-0.4 


4.4 ±0.6 ±0.3 


4.7 ±0.6 ±0.2 


4.6 ±0.4 


5r(5 d -► 7r°7r°) 


< 4.4(90%C.L.) 


< 4.4(90%C.L.) 


< 3.6(90%C.L.) 


< 3.6(90%C.L.) 




< 3.3(90%C.L.) 


< 3.4(90%C.L.) 


< 1.3(90%C.L.) 


< 1.3(90%C.L.) 


Br(B d -► 


< 0.8(90%C.L.) 


< 0.7(90%C.L.) 


< 0.6(90%C.L.) 


< 0.6(90%C.L.) 




< 3.3(90%C.L.) 


< 3.2(90%C.L.) 


< 7.3(90%C.L.) 


< 3.2(90%C.L.) 



For our numerical calculations, we use m,b(mb) = 4.2 GeV for the b quark mass, 
m c (m b ) = 1.3 ± 0.2 GeV, m s (2 GeV) = 110.0 ± 0.25 MeV, (m u + m d )(2GeV) = 9.1 ± 2.1 
MeV. For the decay constants and form factors, we use /„. = 0.131 GeV, f K = 0.160 GeV, 
f B = 0.180 GeV, F B ^ = 0.28, r nK ~ w^r^ = 0.9 [2,18]. For the CKM matrix 
elements, we fix A = 0.2196, and take the central values of | V cb \— 0.0412 ± 0.002, 
| ^ |= 0.087 ± 0.018 [18]. We fix the CP violating phase 7 to 53.6° from global fit- 
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ting of CKM parameters. With the above input parameters, the SM contributions are fixed. 
However, in SUSY the branching ratios depend on several parameters, such as (S^ll, the 
squark mass rriq, and the ratio x = m~/m~. We treat rriq and x as free parameters and con- 
strain (Sjs) ll from the available experimental data. Our main aim will be to obtain limits 
on the complex plane of the flavor changing parameters (5j 3 ) ll for a given x and the squark 
mass. Otherwise mentioned, throughout our analysis we will use 100 GeV for squark mass. 
Constraints on (8jz)LL for other values of rriq can be obtained by rescaling the constraint by 
(rriq/100 GeV) 2 . 

In the theoretical computations of the branching ratios, there are many sources of uncer- 
tainties, mainly arising from different form factors, CKM matrix elements, the annihilation 
contribution and other hadronic parameters of the QCD improved factorization methods. 
We have taken the default values of the above quantities following [2] to get the central 
values of the theoretically calculated branching ratios. However, the results may change if 
one varies any of the above parameters within their allowed range. To take into account this 
fact we conservatively allow an extra 50% variation in the theoretically calculated branching 
ratios from their default values. In the later section we will discuss more on these uncertain- 
ties. The input values chosen here are for an illustration. The main issue of our analysis is 
to show that B — > PP can provide comparable or even better constraints in certain area of 
SUSY parameter space compared with other well known rare processes for B mesons, like 
b — > S7, B — > p7, and the mixing parameters Am^ and Am Bd . 

IV. B -► PP(AS = -1), B S-f AND B° s - B° s MIXING 

In this section we study constraints on the parameter (82%) ll from B — > and compare 
with the constraint from b — > sj. We also comment on how B® — B° s mixing can provide 
information of FCNC parameters in SUSY. Since we are primarily interested in constraining 
the flavor changing parameter (^3) ll, we will present our results on the complex plane of 
the (#23) ll for a given value of rriq and x. 
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FIG. 1. The 2cr allowed region on (623)1,1, complex plane from B — > Ktt with rriq = 100 GeV. 
The solid and dashed bands are the allowed regions corresponding to x = 1 and 4 respectively. 
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-1.5 -1 -0.5 0.5 1 

Re(5 23 ) L L 

FIG. 2. Combined 2a allowed region on (523)ll complex plane obtained from four B — > Kir 
modes as shown in Figure 1. Other notations are same as Figure 1. 

In Figure 1 we show 2a allowed bands on the (<5 23 ) ll complex plane for two representative 
values of x = 1 (solid line) and 4(dashed line) for the squark mass set at 100 GeV. From 
the above Figure with mq=100 GeV, it is clear that the four different decay channels of B 
into K + tt~, K~tc°, K°tc~, K°7r° give almost similar limits and the allowed ranges are almost 
smaller than 1 on the (S23) ll complex plane. There is no significant change in the obtained 
bound even when x varies in the range 0.5 ~ 10. We finally combine the results obtained 
from these four branching ratios and display it in Figure 2. We see that the flavor changing 
parameter (#23) ll is constrained. 

If the squark mass rriq is different than the value 100 GeV used here, one can obtain 
the corresponding constraint on (<5 23 )l L by multiplying a scaling factor m|/(100 GeV) 2 on 
(823) ll- We would like to point out, however, that when making such a rescaling one should 
be cautious about the applicability of the mass insertion method used to obtain the effective 
Hamiltomian. One should only apply the results to regions where (623) ll is not larger than 
of order one. Since the constraint on (S 2 3)ll is already of order one for rrig about 100 GeV, 
for a larger trig B — > PP give weak constraints on the flavor changing parameters. 

We now compare the above constraints with that from b — > sj. We take the current 
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experimental world average result Br(B -> X s -f) = (3.41 ± 0.36) x 10~ 4 [19] from CLEO, 
ALEPH, Belle and Babar. This process is another rare process which is sensitive to new 
physics with AS = — 1 and has been used to constrain new physics beyond the SM. 
The branching ratio of b — > S7 is given by [20,21] 

Br(B - X„) = MB - WB^^ Mmt)!' , (18, 

where Ci2(m&) = cf 2 M (mb) + c^ v '(m&) , g(z) = 1 — 8^ 2 + 8z 6 — z 8 — 24z 4 ln(z), and rj = 
1 — 2/(r, 0, 0)a s (m{,)/37r with /(r, 0,0) = 2.41 [21]. The running of the supersymmetric 
WC's clYu are given by 



cll sy (fi) = r]c™ sv {M susy ), Ci2 sy (/i) = ri 2 cl2 SV (M susy ) + -(r) - ^ 2 )cn SJ/ (M susy ), (19) 
where 77 = (a s (M susy )/a s (m t )) 2/21 (a s (m t )/a s (m b )) 2 / 23 , with 

Cn V (M susy ) = T^i L " ( " ^s(s) - 3M 4 (x)) 

c\r v {M susy ) = - a '< 5 ^ ( - 1|m s (x)) (20) 

The functions M 3!4 (a;) are defined in Ref. [6]. Here we use the leading order result for b — > 57. 
The use of the NLO result [23] will change the LO results by about 10%, which will not 
affect our main conclusion. 

The 2cJ constraints on (823) ll complex plane from b — > S7 for two different values of 
x — 1 and 4 are shown in Figure 3. From this figure we see that as x increases the limit on 
($23) ll from b — > S7 become weaker. 

It was thought that b — > 57 gives the strongest limit on the parameter (823) ll- However, 
in our analysis we find that in certain region of the parameter space the bounds obtained 
from B — > Kn are stronger than the limit obtained from b — > 57. From Figure 4 we see 
that for x = 8 constraints from B — > Kir provide slightly better constraints than that from 
b — > 57. We find that for x less than 8, the constraints obtained from B — > PP are not as 
strong as that from b — > 57. But for x larger than 8, the constraints from B — > PP is more 
stringent. 
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FIG. 3. The 2a allowed bands obtained from b — > 57 on the complex (c>23)ll plane with 
rriq = 100 GeV. Other notations are same as Figure 1. 
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FIG. 4. Combined 2cr allowed region on (623)11 complex plane obtained from B 
band) and b — > 57 (dashed band) for x = 8 with = 100 GeV. 
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One can easily understand the above results from the x dependence of SUSY contribu- 
tions to B — > i^7r and 6 — > 57. From Eq. 7, we find that the SUSY contributions mainly 
appear in parameters a 4 and a 6 through C^J^g. 



susy _ susy , susy _ n LL , ^3 , p f/'TrLZ, , °5 \ 

a — a 4 +a 6 — C 4 + + ^^(^6 "/v" 



(21) 



We have checked numerically that the other contributions in Eq. 7 from SUSY interactions 
are small. 



The parameter a susy linearly depends on ((^Li/m- and depends on x in a complicated 
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but known form which we indicate as a susy = a susy (x). One can easily study how exper- 
imental data may constrain the parameter (t^Li/m- as a function of x by studying the 



behavior of 



R(PP) 



r susy. 



X 



(22) 



a sus y{i) 

The above function shows the behavior of the SUSY contribution with x normalized with 
the value at x — 1. In Figure 5, we show the behavior of R(PP) as a function of x by the 
dark solid line. From the figure we see that R(PP) has a zero value at x ~ 0.37. This 
implies that for x ~ 0.37, no constraint on (623)11 from B — > PP. In the neighborhood 
of this region b — > 57 certainly give a better bound. For x in the range 1 ~ 10, R(PP) 
varies very slowly with x. This almost flat behavior of R(PP) makes the bound on (52z)ll 
to remain nearly unchanged within this range of x. 

5^ , , — , , , — 

4 - 




FIG. 5. R(PP), R(qj) and R(B^ S - B^ s ) as functions of x. The solid, dashed and dotted lines 
correspond to R(PP), R(qj), and R(B^ S — B® s ), respectively. 



While for b — > sj, the SUSY contribution has a different x dependence. In Figure 5 we 
show the variation of 

SUSV { \ 

c 12 y (x) 



R(qi) 



(23) 



cir y (iy 

as a function of x by the dashed line. We see that -R(<n) decreases when x increases. This 
means, for a given m,, the constraint on (^3) ll becomes weaker as x increases. At certain 
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value of x, B — > Kn give better bounds than b — > S7. Numerically we find that when x is 
larger than 8, B — ► i^7r give better constraints on (^23)ll than b — > S7. 

We now comment on the implications of q — q — g interactions on B° s — B° s mixing. The 
effective Hamiltonian for this process is given by [6] 

H e ff = -C m O m , (24) 
a 2 

O m = s a Ll ,b a L s^X (26) 
The functions fs(x) and fe(x) are given by: 



6(1 + 3x) \n(x) + x 3 - 9x 2 - 9x + 17 
6(1 -xf 



M^) = 35 , I 2 ") 



~ 6x(l + x) ln(x) - a; 3 - 9x 2 + 9x + 1 

M*) = ^73^ • ( 28 ) 

After RG evolution we get the WC at m b scale [24]: 

C m (m b ) = rnCmiMwty) , (29) 

^ = / a s (mt) \ 6 / 23 / <Xs ( Msnsy) \ 6/21 Q 

the total AM Bs is given by AM B =2\M? 2 M +M™ sy \ with = Vi(C m /3)M Bs B Bs f Bs . The 

SM contribution M-^ M is given by [25], and M-f^ 27 contains phase from parameter {Sjs) 2 LL - At 
present we have only stringent lower bound on AM Bs which is 14.4 ps -1 at the 95% C.L. [26]. 
Using previously obtained constraints one can obtain the allowed range for AM Bs with SUSY 



contributions. In our numerical calculation we use [18] £ = \JB B J f Bs y B Bd f Bd = 1.16±0.05 
and ^B~B~J Bd = 230 ± 40MeU. 

We find that B — > Kir and b — > 57 data do not constrain AM Bs significantly in a 
large part of the parameter space allowing AM Bs to be much larger than the present lower 
bound. If future experiments can measure AM Bg of the order or less than ~ 48 ps^ 1 , which 
is the limit LHCb can reach, then, one would obtain a much stronger limit on (#23)ll for 
x = 1 and 4 as shown in Figure 6. However we still find that there are regions where 



17 



B — > A"7r and 6 — > 57 give better constraints. To see in which area of parameter space, the 
SUSY contributions to AM Bs become negligible , we study the behavior of relative SUSY 
contribution, R(B® — B®) = C m (x)/C m (l) as a function of x. In Figure 5 we show the 
variation of R(B® — B®) as a function of x by the dotted lines, we see that at x — 2.43 
SUSY contribution to AM Bs becomes zero and no constraint on (623) ll can be obtained. In 
regions around x = 2.43, B — > Kn and b — > 57 can give better constraints even if we assume 
that the future experiment will measure a AM Bs less than 48 ps' 1 . 

1 

.75 
. 5 

j 0.25 

<S 

e 

h -0.25 
-0 . 5 
-0 .75 

-0.75-0.5-0.25 0.25 0.5 0.75 1 

Re (5 23 ) LL 

FIG. 6. The allowed areas (areas inside each curve are excluded if AM Bs is larger than 48 ps _1 ) 
with solid and dashed boundaries corresponding to x = 1 and 4 respectively on the complex plane 
of (<5 2 3)ll- 

V. B -> PP(AS = 0), B p7 AND £ D - 5 D MIXING 

In this section we discuss constraints on SUSY interactions from following processes: 
B — > 7T7T, A" A", _B — > p7 and AM Bd . Constraints obtained from _B — > 7r7r and i? — > XX are 
shown in Figure 7. We find that the constraints on (Si 3 )ll are better than (823) ll obtained 
from B — > A7r decays. The strongest one is obtained from 5~ — > K~K°. 

Although 5 — > p7 has not been measured, the experimental upper bound from BaBar 
Br(B° -> p°7) < 1.4 x 1(T 6 at the 90%C.L. [27] has been used to obtain a constraint on 





\ 


^ 
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FIG. 7. The constraint on (S^ll from B — > 7T7T, B — > if if. For — ► 7r + 7r arc 2cr allowed 
region while for B -► 7r°7r°, B -► K~K° and 5 -► is the limits with 90%C.L. The solid 

and dashed bands (a) and circles ((b)-(d)) are the allowed regions corresponding to x = 1 and 4 
respectively. 
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(8i3)ll- The decay width T(B — > ^7) is given by: 

r(fl - P7) = ^G||^H^||F^(o)| 2 |c 12 (m b )| 2 K + ^) (m| _m ' )3 



m 



(31) 



We use the form factor calculated by T. Huang et al. in Ref. [28], F B ^ P ° = 0.237 ± 0.035. 
The constraints on (S 13 ) LL for x — 1 and 4 are shown in Figure 8. Comparing the constraints 
obtained from B — > nil, KK with B — > p7 we find that similar to the discussions in the 
previous section, when x becomes larger the constraint from B — > K~K° become relatively 
better than 5 — > p7, numerically for the strongest case 5~ — > K~K°, this occurs at re about 
6. 
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FIG. 8. The constraint of (613)11 from B — ► /ry. The allowed area marked by solid and dashed 
line correspond to x = 1 and 4 respectively. 

We now turn to constraint from AM^. This quantity has been experimentally well 
measured with the world average value [18] 0.489 ± 0.008 ps^ 1 . However, the theoretical 
calculation has uncertainty from y/B^f Bd , which is 230 ± 40 MeV. The results are shown 
in Figure 9. It can be seen that AM Bd provides a much stronger constraint on (5i S )ll 
compared with the constraint on (S23) ll from AM Bs . In a large part of the parameter space 
the branching ratios B — > 7T7t, B — > KK provide weaker constraints on (5is)ll unless x is 
close to 2.43 where SUSY contribute to AM# d is zero as can be seen from Figure 5. In the 
region around x = 2.43 B — > K~K° and B — > p7 give better bounds. 
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FIG. 9. The areas marked by "A" and "B" on the complex (S2s)ll plane represent 2a allowed 
regions for x = 1 and 4 respectively obtained from Bd — B<i mixing. 

VI. CONSTRAINTS ON (6 U ) RR 

In the previous discussions, we have concentrated on constraints on (S^ll- We now 
discuss constraints on (o^rr. Due to the different chirality of the new operators the total 
WC's, in the theoretical expressions for B — > PP, are replaced by C|°4^ 6 = Cf^f 5j6 — C^f 5fi , 
where C^ 5 6 have the same form as for C^ 5 6 in Eq. 3 with (5ij) LL replaced by {o^rr. 
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FIG. 10. Combined (B° -► K+tt-, K°7t°) and {B~ -► I^tt , £T°7t-) constraints on (6 23 )rr 
with m 5 = 100 GeV. 
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In Figure 10, we present the combined (B° -> K + n-, K°tv°) and (B~ -> I^tt , fT ?^) 
constraints on (£23)^. The allowed ranges are the same but as a mirror reflection in complex 
plane of (S^rr compared with that of (<%)ll. 
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FIG. 11. b -> 57 and B -+ pj on (5 23 )rr with m,- = 100 GeV. 



In b — * 57 and — > p7 cases, because of different chirality in operators, the Wilson 
coefficient c^ sy is separated from the SM contribution cf^. One needs to replace \ci2{mb)\ 2 
in Eq. 18 and Eq. 31 by |ci2(m b )| 2 = |cf 2 M | 2 + |ci 2 sy | 2 . The allowed ranges are shown in 
Figure 11. We see that they are different than those shown in Figures 3 and 8. It is clear 
that the allowed ranges of (S^rr are constrained to be smaller than 

Bd, s — Bd,s mixing also constrain (S^rr. In this case one just needs to replace ll in 
Eq. 25 by (S^rr and Q m = s a L l^ a L 4lA by Q m = s^^tX- Since (B°\Q m \B°) = 
{B°\Q m \B°), the constraints on (S^rr from Am,B ds are the same as in Figures 6 and 9. 



VII. CONCLUSIONS 



In this paper we studied the constraints on the SUSY flavor changing parameters 
{813) ll,rr and (#23) ll,rr from rare hadronic B — > PP decays. We improved the calcu- 
lations by two folds. Firstly, all the relevant Wilson coefficients are computed at M susy and 
then evolved down to lower scale \x = using Renormalization Group Equation. Secondly, 
we calculated the two body hadronic B decays using QCD improved factorization method. 
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We found that the values of the Wilson coefficients do change significantly under this evolu- 
tion which in turn affect the theoretical computation of B decays. The constraints obtained 
are compared with radiative b — > S7 and B — > p7 decays, and AMs d s . 

We found that each of the different processes discussed can provide interesting constraints 
on the SUSY flavor changing parameters within certain parameter space. For (S23) ll,rr we 
found that as long as x < 7, b — > 57 provides the most stringent limit, however, for x > 8, 
B — > Kit give better constraints. For (5i 3 )ll,rr we observed AM Bd provides a better bound, 
but for x around 2.43, the bound becomes weaker. In this region, B — > K~K° and B — > p7 
can give stronger constraint. 

Before we conclude, we would like to comment on the possible theoretical uncertainties 
arising from different form factors, CKM elements, the hard spectator interactions and also 
from annihilation contributions in QCD improved factorization scheme. As we have already 
mentioned earlier that these hard spectator interactions and annihilation contributions in 
QCD improved factorization method suffer from endpoint singularities and these singular- 
ities have been treated as a phenomenological parameters inducing model dependence and 
additional source of numerical uncertainties. However, the exact estimations of these uncer- 
tainties are not yet known. The uncertainties in the different form factors can change the 
theoretical results by about 40%. We think that the overall uncertainty could be around 
50%. So, in our theoretical calculations we have taken into account this 50% uncertainty 
to obtain the limits. In our calculations we also fixed the SM phase 7 to be 53.6°. When 
SUSY contributions are included, the 7 phase will differ from its SM value. Allowing 7 to 
vary will also change the details of the numerical results. 

However our main objective here was to explore all possible B — > PP decays within 
the context of SUSY and show how in different regions of parameter space the different 
decay modes can constrain the FCNC parameters. The generic feature of our analysis will 
not change significantly by the above mentioned different uncertainties. A coherent study 
of different processes involving 5-mesons carried out here can serve to provide detailed 
information about the SUSY FCNC parameters. We hope in future with the improved 
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theoretical calculations and data will provide a better understanding of SUSY flavor changing 
interactions. 
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